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Abstract 

We construct examples of compact and one-ended constant mean curvature surfaces with 
, large mean curvature in Riemannian manifolds with axial symmetry by gluing together small 

spheres positioned end-to-end along a geodesic. Such surfaces cannot exist in Euclidean space, 
but we show that the gradient of the ambient scalar curvature acts as a 'friction term' which 
permits the usual analytic gluing construction to be carried out. 



1 Introduction 

> 

^T) I Background. The study of constant mean curvature (hereafter CMC) surfaces in R'^, or more 



generally in three-dimensional Riemannian manifolds, is a well established field of Riemannian 

. geometry and the literature concerning the construction and properties of such surfaces is enormous. 

. One particular method for constructing CMC surfaces is by analytic gluing techniques. These go 

■ back to the work of Kapouleas [6],[7j and have been further developed by many others, including 

the first author with Pacard [T] , [2] and the second author with Pacard and also with Pollack [S] , [TT] . 

X ■ (See [12] and [l6] for surveys about the current state of this approach.) The general idea here is 

' 

' to take connected sums of simple surfaces (for instance the classical examples of nonminimal CMC 

surfaces in — the sphere, the cylinder, and the one-parameter family of Delaunay surfaces 
of revolution) to produce more general CMC surfaces of finite topology, both compact and non- 
compact. Such constructions work by finding a surface whose mean curvature is nearly CMC 
embedding of the connected sum and perturbing it to have constant mean curvature. 

There are numerous constraints on the structure of CMC surfaces of finite topology in M.^. For 
example, Meeks [13] proved that any end of a complete Alexandrov-embedded CMC surface in 
is cylindrically bounded; while Korevaar, Kusner and Solomon [8] improved this by showing 
that any such end converges exponentially to one end of a Delaunay surface. Furthermore, the 
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possible directions of the axes of these ends are also subject to limitations, as well as the flexibility 
to change these directions within the moduli space of all such surfaces, see [1],I5] and [3]. These 
limitations are phrased in terms of a certain flux that was discovered by Kusner. The flux is a 
vector associated to any closed loop in a CMC surface; it is constant under deformations of this 
loop, and in fact only depends on the homology class of this loop in the surface. There is a flux 
associated to a simple positively oriented loop around each asymptotically Delaunay end, which 
depends only the direction of the axis and the neck-size of the limiting Delaunay surface. The 
homological invariance also shows that the sum over all ends of these limiting fluxes must vanish, 
which is a global balancing condition for the entire CMC surface. One immediate consequence is 
the non-existence of a complete Alexandrov-embedded CMC surface in with only one end. The 
flux also provides useful local information. For instance, it is a crucial ingredient in the gluing 
constructions mentioned above, since as we explain more carefully below, one must choose the 
initial approximate CMC configurations so that the fluxes are almost constant across the necks of 
the connected sums. 

We turn now to a newer theme in this subject, namely the study of sequences of CMC sur- 
faces with mean curvature tending to infinity in an arbitrary 3-manifold. Examples in include 
sequences of dilations of the complete CMC surfaces with k asymptotically Delaunay ends. Such 
sequences of surfaces 'condense' onto one-dimensional sets, here a union of half-lines meeting at a 
point. This seems to be a general phenomenon: Rosenberg [18] has shown that if S is a closed CMC 
surface in an arbitrary (compact) 3-manifold Af, with sufficiently large mean curvature H, then 
M \ T, has two components, and the inradius at any point in one of these components is bounded 
above by C/H. In other words, S looks like a tube around some (presumably 1-dimensional) set 
7. If Tij is any sequence of CMC surfaces with mean curvature Hj — > 00 condensing to a curve 
7, then a formal calculation (assuming that the supremum of the pointwise norm of the second 
fundamental form is of the same order as Hj) shows that 7 is a geodesic, or at least a union of 
geodesic arcs. This leads to the following central question. 

Question: What are the possible condensation sets 7 in a 3-manifold M for sequences 
of CMC surfaces Sj with mean curvatures Hj y 00. 

The obvious guess is that a condensation set is some sort of network of geodesies. Based on the 
examples of dilated CMC surfaces in M'^, one expects each edge of this geodesic network to have 
a 'weighting' which carries information about the Delaunay parameters of the CMC tubular piece 
which converges to that edge. This is far from being proved, but there are some very partial results. 
Motivating this conjecture and as a first step on it, the second author and Pacard [lOj proved that 
if 7 is any closed geodesic in M which is non-degenerate (in the sense that its Jacobi operator is 
invertible), then most geodesic tubes of sufficiently small radius about 7 can be perturbed to CMC 
surfaces with large H. The present paper undertakes a next step toward this conjecture. We prove 
that for very special (non-compact) ambient 3-manifolds, there do exist sequences of one-ended 
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CMC surfaces condensing to geodesic rays and sequences of compact CMC surfaces condensing to 
geodesic intervals. However, for reasons we explain below, it is not clear whether the fact that 
these limiting curves are geodesies is their most relevant feature. Our construction works in higher 
dimensions too, i.e. we construct sequences of CMC hypersurfaces condensing to a ray or interval. 

There is an important feature of this condensation question which has not been mentioned 
yet. The first result about CMC surfaces of high mean curvature was due to Ye [20J in the early 
1990's. He proved that if p is a non-degenerate critical point of the scalar curvature S of M"^^, 
then geodesic spheres around p with small radius may be perturbed to CMC surfaces with large 
H. He also showed the converse: assuming bounded eccentricity, sequences of CMC spheres with 
H ^ oo converge to a point p where \7S{p) = 0. A more recent paper by Pacard and Xu [14J 
considers the same problem in manifolds with constant scalar curvature and proves that there is 
a secondary curvature function whose critical points regulate the location of these small CMC 
spheres. Therefore the role of the scalar curvature of the ambient manifold in the question of CMC 
surfaces condensing to one-dimensional sets must be addressed. 

The CMC surfaces constructed in this paper are perturbations of collections of small spheres 
joined together by even smaller catenoidal necks, all arranged along a curve 7. The sizes of the 
catenoidal necks are quite small compared to the radii of the spheres, but quite strikingly, these 
neck-sizes must vary along this chain of spheres. Indeed, it is precisely because these neck-sizes 
decrease that the surface eventually 'caps off' to an end rather than continuing. This phenomenon is 
caused by a flux formula that involves the gradient of the scalar curvature of the ambient manifold. 
Indeed, unlike in Euclidean space, the difference of the fluxes computed on two loops which are close 
to one another need not vanish, but may be computed in terms of a surface integral involving VS. 
(This is completely analogous to the generalized Pohozaev identity discovered by Schoen which 
arises in his construction of metrics of constant positive scalar curvature [ISj.) In our setting, this 
shows that the difference between successive neck-sizes can be expressed in terms of the gradient 
of the scalar curvature along the axis connecting these two necks. Thus, in some sense, V/S acts as 
a friction term. By contrast, the scalar curvature of M does not play a role in the location of the 
1-dimensional condensation set in [10]; the only important feature there is that the condensation set 
is a closed non-degenerate geodesic. This is almost surely because the CMC surfaces constructed 
there are nearly cylindrical. 

We now describe this more carefully. Let S be a hypersurface with constant mean curvature H 
in (M"+^, (?). Suppose that U and W are open sets in S and M, respectively, such that dW =U\JQ 
for some hypersurface-with-boundary Q. If there happens to exist a Killing fleld V on M, then the 
flrst variation formula for the area of li with the volume of W flxed relative to the one-parameter 
family of diffeomorphisms for V gives that 



here u is the unit normal vector field of V in E and N is the unit normal vector field of Q in M. 




(1) 
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The flux itself is defined as 

/ g{u,V)-H [ g{N,V) (2) 

where 7 is a curve in S and Q is any surface in M with dQ = 7; this integral is independent of the 
choice of Q. Thus ([T]) is the statement that the flux depends only on the homology class of 7 in S. 

As already noted, these flux integrals determine when an approximately CMC surface can be 
perturbed to be exactly CMC. This can be explained more concretely as follows. Suppose flrst 
that M = W"^^ and let S consist of a collection of spheres of radius r (hence mean curvature 
n/r) connected to each other by small catenoidal necks. Let lA be one of these spheres with two 
small spherical caps removed where the necks are attached, Q the union of two disks capping these 
boundaries and W the slightly truncated ball enclosed by ^ U Q. Then ([1]) becomes 

/ g{u,V)-- [ g{N,V)= V r^~%V, + 0{r^-'e^), (3) 

■^9" all necks 

where Vi is the unit vector pointing from the center of sphere in question to the i*^ neck, rsi is the 
width of this neck and e := maxjjej}. If S were exactly CMC, the left hand side would necessarily 
vanish. If S is not exactly CMC, then it is a fundamental fact that in order to find a nearby CMC 
surface, it suffices that the leading term on the right hand side of ([3]) must vanish for each spherical 
region U in S. If this condition is satisfied, the approximate CMC surface S is called balanced. 
Note that it is impossible to have a balanced approximately CMC surface where some sphere has 
only one spherical neighbour. 

When the ambient manifold (M"+^,(7) is arbitrary, one can form approximate CMC surfaces 
with large H by attaching together some large collection of geodesic spheres of very small radius r. 
There are (usually) no Killing fields, but we can use the approximate Killing fields corresponding to 
translations and rotations in Riemann normal coordinates based at the center of any one of these 
spheres. Formula ([T|) now becomes 

/ g{v, V)-H f g{N, V) = V] r^'-^SiVi - Cr''+'^VS{p) + 0{r''-^e^) + 0(r"+^) (4) 

Jq ,u necks 

where VS* is the gradient of the scalar curvature of M and C is some explicit dimensional constant. 
Note that when applied to a sphere with no neighbours, hence with all = by default, this gives 
Ye's condition that the right hand side of (HD must vanish like r""'"'^. 

The main point in this paper is to exploit the contribution of VS* in Spheres of radius r 
are joined by necks of width re, where by this same formula it is natural to assume that e = 0{r^); 
these configurations are arranged in such a way that the leading term on the right in ^ vanishes. 
The perturbation argument producing a nearby CMC surface is then not so different than the one 
in Euclidean space. 

Description of the surfaces. Two specific examples of CMC surfaces in M exhibiting markedly 
different properties from those occuring in Euclidean space will be produced in this paper. We 
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shall make extremely strong assumptions about the geometry of {M,g) in order to simplify the 
calculations, which even so are still quite lengthy. Thus our result should be regarded as a model 
for what should happen in more general manifolds, though that would take considerably more 
work. Our surfaces will consist of geodesic spheres of small radius r arranged along a geodesic 
segment or ray 7 C M and joined by suitably scaled pieces of catenoids. In the first example, 
some large number of spheres, on the order of 1/r, are glued together so that the resulting surfaces 
are embedded and compact; the second is embedded and complete, and is built from some large 
number of spheres joined at one end to a half Delaunay surface of small neck size. In either case, 
there is a terminal spherical component which has only one spherical neighbour. 

We make the following assumptions about {M,g). First, let 7 be a geodesic segment or ray in 
M and assume that there is a neighbourhood of 7 in which the metric g is axially symmetric, i.e. 
invariant with respect to rotations about the axis 7. Thus Fermi coordinates around 7 identify this 
tubular neighbourhood with [0,L] x D (where L = 00 is allowed and D is & disk in M"^^), and 

g = dt^ + A{t)5- (5) 

here 5 is the standard Euclidean metric on D and A{t) is a smooth, strictly positive function of 
the arclength t along 7. The scalar curvature of 5 is 5 := A~'^{—nAA + A^). Further 

assumptions on A depend on whether we wish a finite-length or one-ended CMC submanifold. 

1. In order to construct a finite-length surface, assume that A is an even function of t so that the 
reflection t ^ —t induces an isometry of the tubular neighbourhood of 7. Assume furthermore 
that f = is a non-degenerate local maximum of S. 

2. In order to construct a non-compact one-ended surface, assume that when t > 0, the scalar 
curvature is negative and increases monotonely to 0, and that \S{t)\ < Ce"* for some a < 0. 
(One function which satisfies this is A{t) := 1 -|- e~*.) 

These assumptions significantly reduce the complexity of the perturbation argument, but allow 
for the one feature which allows for this new behaviour of CMC surfaces, namely that VS* points 
along 7. By our assumptions, however, the geodesic 7 is also an integral curve for VS, and it 
is unclear which of these two features is the crucial one. One basic question we leave open is 
the geometric characterization of these condensation curves in more general ambient geometries. 
We expect new and interesting behaviour to occur when VS" no longer is required to point along 
geodesies. 

Our main result can be expressed as follows. 

Theorem. Let M he a Riemannian manifold with the special features described above. 

• Let L be a finite segment of the geodesic 7, where the parametrization of 7 is such that the 
reflection t ^ —t is an isometry in some neighbourhood of L. Then there exists an r^ > 
so that for every < r < ro, there is a CMC surface which is a small perturbation of a 
surface constructed by gluing together length{L)/r spheres of radius r with centers lying on 7. 



5 



• Let I be a ray of the geodesic 7. Then there exists an r^ > so that for every < r < r^, 
there is a CMC surface T,^^ which is a small perturbation of a surface constructed by gluing 
together a number 0(l/r) spheres of radius r with centers lying on 7, together with an end 
of a Delaunay surface whose axis lies along 7. 

The rescalings of and by the factor 1/r converge as r ^ to an infinite or semi-infinite 
string of spheres of radius 1 with centers arranged along a segment or ray in M". The precise mode 
of convergence wiU become clear in the course of the proof. 

As already noted, the proof roughly follows the proofs of analogous gluing theorems for CMC 
surfaces in Euclidean space: an approximately CMC surface S is deformed via small normal defor- 
mations which are parametrized by functions on S, which transforms the problem to one of finding 
a solution of the constant mean curvature PDE. One difficulty is the fact that one expects to find 
a solution only when all parameters are very small, but this means that the geometry of S, and 
hence the PDE which must be solved, are very degenerate. In addition, the Jacobi operator (i.e. 
the linearized mean curvature operator) on S has small eigenvalues generated by the nullspaces 
of the Jacobi operators on each spherical and catenoidal component. The PDE is first solved on 
the finite codimensional orthogonal complement of this approximate nullspace. By repositioning 
the various components of this approximate solution one can show that it is possible to kill the 
remaining finite dimensional piece too, provided the map carrying the 'repositioning parameters' 
to the right hand side of the flux formula dH) is an isomorphism. One must also keep careful track 
of the dependence on the neck-sizes e and radii r in all of this to guarantee that the estimates 
controlling the existence of the solution of the constant mean curvature PDE are strong enough. 

Acknowledgement. The authors wish to thank Frank Pacard for interesting discussions during 
the course of this work. 

2 Preliminary Geometric Calculations 

The approximate solutions constructed in Section [3] are assembled from small geodesic spheres 
centered on points of the geodesic 7 in M connected to one another by small catenoidal necks. It 
is most convenient to use geodesic normal coordinates centered at points of 7. Since the ambient 
metric 5 is a second order perturbation of the Euclidean metric in these coordinates, the first step 
in every estimate is to perform the computations for a Euclidean metric; the second step is to 
incorporate the perturbations coming from the metric into the estimates. In this section we derive 
various expansions of the mean curvature and other geometric quantities. 
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2.1 Geometry of Surfaces in a Geodesic Normal Coordinate Chart 

If p is any point in a Riemannian manifold (M, g) , then in terms of geodesic normal coordinates 
centered at p, 

g:=g + P ■= {Sij + Pijix)) dx' (g) dx^ 
where g is the Euclidean metric and P is the perturbation term. It is well known that 

Piji^) = IJ2 Rajmrnx^x"" + RiljmA^)^'^"^'' + 0{\\xf). (6) 

l,m l,m,n 

where Rijki ■= Rm(^, ^) and Rijki-m ■= V_|__Rm(^, ^) are components of 

the ambient Riemann curvature tensor and its covariant derivative (the ambient covariant derivative 
is denoted V). 

Suppose that S is a surface in M. The following results provide expansions for various geometric 
quantities of S in terms of P. Here and in the rest of the paper, let h, T, V, A, N, B, H be the induced 
metric, Christoffel symbols, covariant derivative, Laplacian, unit normal vector, second fundamental 
form, and mean curvature of S with respect to the metric g, and let h,T,V , A, N, B, H be these 
same objects with respect to the Euclidean metric. Near a point x G S, let {Ei,E2} be a local 
frame for TS induced by some coordinate system and denote by Y := ^"'^ the position vector. 
Define 



Poo 


:= P{N,N) 






:=P{N,Ej) 






= PiE,,E,) 






:= \{{E,P){E,,Et) + {EjP){E,,Et) 


-{EtP){E,,E,)) 


Vijo 


■■=l{iEiP){Ej,N) + {EjP){Ei,N) 


-{NP){Ei,Ej)). 



Straightforward geometric calculations now yield the following results. 

Lemma 1. The induced metric oJY, and the associated Christoffel symbols are given by 

hij = hij -\- T^ij and ^ijk ~ ^ijk ~t~ ^ijk ~t~ ^OkRij • 

The normal vector of E satisfies 

{l + VoQ-h^VoiVojf^ ' 

The second fundamental form of S satisfies 

= (1 + Poo - h'^VoiVoj) + Pvo - h^^PokV.,i 

{l + Voo-h^^VoiVoj)'^' 
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Lemma 2. The induced metric ofE satisfies 

The normal vector of S satisfies 

N={1- \Voo)N - h'^VojEi + 0{\\Y\\^) . 

The second fundamental form and mean curvature of S satisfy 

Bij = (1 + ^Voo) Bij + Vijo + Bij{Y, B, N, ^1,^2) 
H={1 + \Voo)H + U^VijQ - B'^Vij + H{Y, B, N, Ei,E2) . 

where Bij and H are functions satisfying 

max\Bij{Y,B,N,Ei,E2)\ + \HiY,B,N,Ei,E2)\ < C\\Yf{l + ||y||||i||) 

for a constant C depending only on the curvature tensor of the ambient manifold at the center of 
the normal coordinate chart under consideration. 



The mean curvature of S may now be computed as 

H={1 + iRm(iV, Y, N, Y) + VyRm(7V, Y, N, Y))H 

- {\-Rm{Eu y, Ej,Y) + |VyRm(£;i, F, Ej,Y))&^ 



(7) 



- |Ric(y, A^) - |VyRic(r, AT) + j^V^Ric(y,y) - iV^Rm(Ar,y,Ar,y) 

+ 0{m\^ + \\Y\\\\B\\))- 
The third hne here contains the largest terms coming from the ambient curvature. 

2.2 Mean Curvature Calculations in Euclidean Space 

Let S be a surface in Euchdean space. Choose a function / : S ^ M and define fif : T, ^ M""*"^ 
to be the normal deformation of S by f{p). The mean curvature operator / 1— > iJ[/i_^(E)] with 
respect to the Euclidean metric decomposes as 

H[fif{^)]=H + l{f) + Q{f) (8) 

where JC.{f) := A/+ ||-B|p/ is the linearized mean curvature operator with respect to the Euclidean 
metric and is the quadratic and higher remainder term. The second fundamental form can be 
similarly expanded as B[fif{^)] =B + M^\f) + B^'^\f). We now derive expansions for 
B^'^\f) and Q{f) in terms of /. Although these results are fairly standard, it is important is to 
track the dependence on in the various error terms appearing in the expansions. 

We first expand Bf := B[fif{'£)] and Hf := H[fif{Il)] in terms of / and extract the constant, 
linear and higher-order parts. Introduce 

I3st :=hst- fBst P'' := [inverse of /3] D := {l + (3"" (3^' f,kf,ihijf'^ 
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where a comma denotes ordinary differentiation in tlie coordinate directions. After some work, one 
finds that the Euchdean induced metric hj, its inverse, and the Euclidean unit normal vector Nj 
of /^/(S) are 



L>2 



The second fundamental form Bj := B[flf(Y,)] can thus be expressed in terms of / as 

[^f]ij = { ^ij + /;u ~ fBiBkj + 13^^ f,i{f,iBjk + fjBik + /Bj^.i) J (9) 



where a semicolon denotes the covariant derivative of S with respect to hst ■ We expand the inverse 
of the induced metric as 

hj := °0 + 2fB'^ - ■^f3'"'P^''P'Pp"^Uf,ihmphng + r?'"V"/i™n 

where the remainder r/*-' in := h^^ + /i?*-' + rf^ satisfies \rf^\ = 0{\f\^\\B\\^). Now taking the 
trace of ([9]) with respect to nj yields the mean curvature Hf := H{p,f{T,) which is 

DxHf = H + Af+ \\Bff + 2fB'if,i, - 2fTts{B^) + h'^ P^' f,i{V,iBjk + fB,k;i) 

+ (^-^P'"'P'''P'Pp'''f,kf,ihmphng + r?*"r?^" (Bij + f,,, - fBfBkj) (iq) 



All of this is summarized in the following lemma. 
Lemma 3. The linear parts of Bf and Hj are 

[B^'Hmr-=Uj-fBfBkj 

£(/) ■.= Af + \\Bff. 
The quadratic remainder parts of Bf andHj are 



[B^^Kf)]ij ■■= ^{fAk + fjB^k + fBjk;i) + - l) i^^J + - fB'Bkj) 
2fB''f-ij - 2f^TTs{B^) + W^P^'f,i{2f,iB,k + fB.u.i) 
+ (-^l^'""P'''P^''P''f,kf,ihmphn, + V'"'r^'''h^n) {B^j + f.,ij - fBfBkj) 



Qif) ■= ^ 



+ (2fB'^ - ^P'"'f3^''P''Pf3'''f,kf,lhmphnq + r^'^^rf^hmr^ f3^'f,i{2f^,Bjk + fBjk;i 



+ {H + Af + \\Bff) 
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The quadratic parts of both Bf and Hf are unwieldy, but only basic structural facts about them 
are needed in the sequel. To simplify matters, we suppose that |/|||-B|| + ||V/|| <^ 1, which will be 
justified later on. Furthermore, both B^'^\f) and Q{f) can be expanded into a sum of terms which 
are each hnear combinations of the coefficients of the tensor /* • (V/)'^-' (V^/)'^'^ iJ'^'ig) (V-B)*^™, 
where i,j,k,l,m are positive integers such that k + m < 1 and i + 1 = k + l + 2m (i.e. the number of 
times the function / appears is smaller by one than the sum of the number of covariant derivatives 
and the number of occurrences of the second fundamental form) . Consequently the dominant terms 
in B^'^\f) and Q{f) are 0(1) linear combinations of components of 

{Vff^B, fVf^B\ fVf^VB, fB\ fV^f^B and {Vff^V^f. 

The following estimates are now straightforward consequences of this discussion. 

Lemma 4. Assuming that ||i?|| + ||V/j|| ^ 1 for i = 1,2, the quadratic remainders in the second 
fundamental form and mean curvature satisfy 

\B(^)if,)-B(^\f,)\ + \Qif,)-Qif,)\ 

< C\fi - /2I • max {\fMBf + \\Vfi\\\\Bf + ||V/i || || Vi?|| + \\V^ fi\\\\B\\) 
+ C\\Vfi - V/2II •maxd/.lpf + ||V/,||p|| + \fi\\\VB\\ + \\V^fi\\) 
+ Clival - V2/2II • max + ||V/,||) 

+ Clival - V2/2II • max \\Vfif + C||V/i - V/2II • max || V/,|| ||VVi|| 

i i 

where C is independent o//i,/2 and \\B\\- 

2.3 Mean Curvature Calculations for a Perturbed Background Metric 

Consider now a surface := jJifiTi) deformed by the amount / in the direction of the Euclidean 
normal to S. Working again in a geodesic normal coordinate system centered at some point of 7, 
we now decompose the mean curvature operator / ^ Hf := ^^[5^/] as 

Hf = H + C{f) + Q{f) + -H{f) (11) 

into constant, linear and quadratic remainder parts plus a small error term. 

The key is to substitute the tangent vector fields [Ef]i := y^j3ijEk + /,jiV of /i/(S) (here Ei,E2 
are tangent vectors for S), the Euclidean normal vector field Nf := j){N — (3^^ f^^^Ej^ of fifiTi), 
and the position vector field Yj := y + fN of /U/(S) relative to the center of the normal coordinate 
chart, as well as the expressions for H and B into the formute from Lemma [2j This yields 

H={1 + RiU))H + [R2{f)hB'' + R{f) + nf) , (12) 

where Ri, R2 and R are first-order differential operators and Ti.{f) := 7i{Yf,Bf,Nf, [Ef]i, [Ef]2). 
As before, the precise structure of these operators is not important, though we still must estimate 
their dependence on /, y and B. 
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First, by examining the expansions for Rg in terms of f,Y and V/, and for D and (3 in terms of 
fB and V/, one deduces that Rs has an expansion into constant, hnear and quadratic remainder 
terms of the form Rs{f) := R^^ + Ri^\f) + Ri^\f) where 

Rf\f) ■■= lRmiN,Y,N,Y) 
R\'\f) := -iRm(iV,F,i?,-,y)^^V,fe 
[RfHf)]ij:='eRHEi,Y,Ej,Y) ^^^^ 
Ri'^ (/) := i {Rm{Ei, Y, Ej,N) - Rm{Ek, Y, E^, Y)B^) f + iRm(iV, Y, E^, Y)fj 
+ {lVYRrn{Ei,Y,Ej,N) + j^RTa{Ei,Y, Ej,Y)) f 

(2) 

for s = 1,2. Here R],q is a sum of quadratic and higher expressions in the components of 
jfcy®{n-fc)(^(V/)®' for various n > 2, k < n and / > whose coefficients are bounded by curvature 

(2) 

quantities, while R). [ can be expanded to any order in a power scries in the components of fB and 
V/ whose coefficients are bounded by curvature quantities. One finds a similar expansion for -R(/) 
into constant, linear and quadratic remainder terms of the form R{^f) := R^^^'^ + R^^\f) + R^'^\f) 
where 

:= -iRic(y,iV) - iVyRic(y,iV) + ^V^RiciY,Y) - ^Rm{N,Y,N,Y) 
R^^\f) ■■= -(iRic(iV,iV) + iVyRic(iV, AT) + ^Vj^Ric{Y,N) - |V^Ric(iV, Y))/ 

+ (iRic(y, Ei) + i VyRic(y, Ei) - ^ Vfi,Ric(y, y)) h'^fj (u) 

+ {lVj;,Rm{Ei, y, TV, Y) + ^ V£;,Rm(iV, Y, N, Y)) h'^ fj 
(/) := V/, y) • {fB, <Jf) 

- (2) . 

for s = 1,2. Here Rq is a sum of quadratic and higher expressions in the components of 

various n > 1, k < n and I > whose coefficients are bounded by 
curvature quantities, while R^ can be expanded to any order in a power series in the components 
of fB and V/ whose coefficients are bounded by curvature quantities. Therefore the following 
estimates hold. 

Lemma 5. The constant and linear parts of Rs{f) for s = 1,2 and R{f) satisfy 

\RfHf)\<C\\Yf 

\R^\f)\<C\\Y\\{\f\ + \\Y\\\\Vf\\) 
\R(^\f)\<C\\Y\\ 
\B'\f)\<C{\f\ + \\Y\\\\Vf\\). 

The quadratic and higher parts of Rs{f) for s = 1,2 and R{f) satisfy 

\Ri'Hfi) - Ri'\f2)\ < - /2I + ||y||||V/i - V/2II) • max (|/,| + ||y||||V/.||) 
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\R(^){f,) - R('\h)\ < C\h - /2I • max ||V/i|| + C||V/i - V/2II • max \fi\ 

+ C{\h - /2I + nil|V/i - V/2II) • max (|/,| + ||y||||V/,||) . 

In these estimates, C is a constant depending only on the curvature tensor of the ambient manifold 
at the center of the normal coordinate chart under consideration. 

One now substitutes the expansions for R{f) and R{f) along with the expansions for Hf and 
Bf in terms of / into equation (jl2jl and extracts the various parts. That is, by performing these 
substitutions, one finds 

H:= (l+i?f))F + i?W+<).S(°) 

CU) := (1 + Rf^)CU) + R^'Hf) + R^'^ ■ B^'Hf) + R^'Hf) ■ B (15) 
Qif) := (1 + Rif))Qif) + R^^Hf) + • + ii« (/)£(/) + Rif) ■ B^'Hf) 

where [R{f)]ij '■= Ri{f)hij + [i?2(/)]ij- Moreover, the following estimates hold. 

Lemma 6. The operator C satisfies 

\c{u) - c{u)\ < c{i + ||y||||^||)(|u| + ||y||||Vn|| + ||yf nv^nH) (le) 

where C is a constant depending only on the curvature tensor of the ambient manifold at the center 
of the normal coordinate chart under consideration. 

Lemma 7. Under the assumption that |/j|||i3|| + ||V/j|| <C 1 for i = 1,2, the quadratic remainder 
Q{fi) satisfies 

|Q(/i) - Q(/2)| < C|/i - /2I • max {\fi\\\Bf + ||V/,||||^f + ||V/,|| ||VS|| + HVVdlPII) 
+ CIIV/1 - V/2II • max(|/i|||5||2 + ||V/i||||^|| + \fi\\\VB\\ + fi\\) 
+ Clival - V2/2II • max {m\\B\\ + \\<J fi\\) 

+ Clival - V2/2II • max ||V/,f + C7||V/i - V/2II • max ||V/i|| HV^/ill 

i i 

+ C\fi - /2I • max WVfiW + C7||V/i - V/2II • max 

i i 

+ - /2I + \\Y\\ llV/i - V/2II) • max {\fi\ + \\Y\\\\V f,]\) . 

In these estimates, C is a constant depending only on the curvature tensor of the ambient manifold 
at the center of the normal coordinate chart under consideration. 

Finally, further straightforward calculation leads to the remaining estimate for TC. 

Lemma 8. Under the assumption that |/|||i?|| + ||V/|| <^ 1 then the error term Tiif) satisfies 

Mfi) - n{f2)\ < c\\Y fii + ||y||p||)(|/i - /2I + ||y||||v/i - V/2II + ||yf liv^/i - V2/2II) 

where C is a constant depending only on the curvature tensor of the ambient manifold at the center 
of the normal coordinate chart under consideration. 
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3 The Approximate Solutions 



Wc now construct two families of approximate CMC smfaces. The first family consists of finite- 
length surfaces invariant under the reflection t t— > — t constructed by gluing together K small 
geodesic spheres of radius r along the (t, 0, 0) geodesic with small interpolating necks. Here K 
is approximately 1/r so that the surface fills out a region along 7 of bounded length which does 
not tend to with r. The second family consists of one-ended surfaces in an asymptotically flat 
Riemannian manifold. These are constructed by taking a configuration of K spheres as above and 
then attaching a semi-infinite Delaunay surface to the last sphere. These two families are denoted 
S^((7, 5) and S^^((7, 5), respectively. These depend on parameters (Ji, cr2, . . . and (5i, (52, . . . which 
govern the precise locations of the component spheres and necks. For brevity, we often just write 
Sr(cr, (5) for either family when the context is clear or does not matter. 

Remark: Starting from this point, our presentation will be phrased in terms of two-dimensional 
surfaces S contained in a three-dimensional Riemannian manifold M. This is done for the purpose 
of simplicity; however, everything that follows can be easily adapted to the (n -|- l)-dimensional 
setting. 

3.1 The Finite-Length Surface 

Let 7 be the (i, 0, 0) geodesic and (with slight abuse of notation) also the arc-length parametrization 
of this geodesic given by 7(t) := (t, 0, 0). Introduce the small radius r and the (much smaller) 
separation parameters Uk- Let = and tk '■= 2kr + Yl'i=oi'^l) ^^'^ P±k '■= (i^fejOjO). 
The geodesic spheres that will be used in the gluing construction are S±k ■= dBr{p±k)- Also let 
p\: ■~ ^i^k + r + crfe/2) be the point half-way between Sk and Sfc+i and let := ^{tk ± r) be the 
points in Sjfc n 7 . Define and p^_j^ in a symmetrical manner. In the definitions above, the index 
k ranges from zero to K. 

The construction of the first surface consists of three steps. The first step is to replace each 
with Sk which is obtained from Sk \ {p^ or Sq by a small normal perturbation designed to 
make Sk look more like a catenoid near p^. The next step is to find the truncated and rescaled 
catenoids that fit optimally into the space between Sk and its neighbours, the precise location of 
which is governed by displacement parameters Sk- The final step is to use cut-off functions to 
glue each Sk smoothly to its neighbouring necks. The result of this process will be a family of 
surfaces that depends on r and the parameters cri, (T2, • • • and 81,62, ■■ ■■ A number of additional 
small parameters and will be introduced below and it will be shown how these depend on 
the a and 5. Denote by e := maxfcjefe, e^^, ej^} below and in the rest of the paper. 

Step 1. Let jCs '■= ^ + ll-Bp be the linearized mean curvature operator of Sk with respect to 
the Euclidean metric and let Jk ■ Sk ^ ^ he the smooth function in the kernel of £5 that is 
cylindrically symmetric with respect to the axis defined by the geodesic 7 and normalized to have 
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unit L^-norm. (It is defined by taking the correct multiple of the Euclidean normal component of 
the translation vector field ^.) 

Now introduce small positive scale parameters that have yet to be determined and define 
Gk '■ Sk \ {p^} for k = 0, . . . , K — 1 as the unique solution of the equation 

Cs{Gk) = e+6+ + e-6. + AkJk (17a) 

where 6± denotes the Dirac 5-function at and the real number Ak is chosen to ensure that the 
right hand side of (jl7ap is L^-orthogonal to Jk- Also, let G'q be the unique solution of the equation 

Cs{G'k) = e^6- + A'^.Jk (17b) 

where A'j^ is chosen to ensure that the right hand side of ()17bp is L^-orthogonal to Jk- Note that 
l^fcl ^ Ge for A; = 1, . . . , and = by symmetry. 

To complete this step, introduce another small radius parameter yet to be determined, 
and define Sk as the Euclidean normal graph over Sk \ [-B^j./2(Pfc ) U B^^/2iPk)] ^^^^ generated 
by the function rGk- Also define the terminal sphere S'j^ as the Euclidean normal graph over 
Sk\B^.j^I2{p^) that is generated by the function rG'^ as well as its symmetrical counterpart under 
the t ^ —t symmetry. 

Step 2. Coordinatize a neighbourhood olp\ using geodesic normal coordinates centered at p\ and 
scaled by a factor of r. Let the coordinate map be V'fc : Br{p\) C {—R',R')xBri{Q) CMxM^, 

where R,R' are appropriate radii (one should think of i? = 0{r) and R' = 0(1)). Note that the 
R coordinate corresponds to a translation of the scaled arc-length coordinate along 7 and 7 itself 
maps to the curve ^ (x'^, 0) with mapping to the origin. The gluing procedure that will now 
be described applies to any pair of perturbed spheres Sk-, Sk+i-, including the last pair Sk-i,Sk- 
The images of 5*^+1 and Sk under the coordinate map, at least near the origin, can be represented 
as graphs over the M? factor of the form {{x^,x) : x^ = F^fiix)} where + and — correspond to 
ipk{Sk+i) and ipkiSk) respectively. One can check that the Taylor series expansions for Gk near 
p^ and for ipk imply corresponding expansions for F^f^ near x = 0, which results in the fact that 
tpkiSk-^^i) is the set of points 

^° = ^.m(^) 5 + ^fc+ifei + ^k+i log(lkll)) + (^(Ikf ) + Oie\\xf) 
while ipkiSk+i) is the set of points 

= FU"^) ■■=-Yr- '^^"^ + ^^S{\m + 0{\\xf) + 0{e\\xf) . 
Here, and are constants. 

The interpolation between ipkiSk+i) and ipki^k), will be done using a standard catenoid that 
has been scaled by a factor of Sk and translated by a small amount along its axis. The x'^ > end 
of such catenoid is given by 

= Fneck(^k,dk;x) := arccosh(||x||/efc) +efc4 
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= log(2) - log(£ifc)) + Ek log(||x||) + £fc4 + C(4/||xf ) 

near the origin, where dk is the translation parameter. The x° < end is given by 

= Kecki^k, dk; x) := - £fc arccosh(||a;||/efc) + eudk 

= - £fc(log(2) - log(efc)) - 8k log(||x||) + £kdk + Oiel/\\xf) 

near the origin. Optimal matching of these asymptotic expansions of the catenoid with those of 
tpkiSk+i) and ipkiSk) given above then requires 

and 



,+ 



ak = Akisk) := rsk (^2(log(2) - log{ek)) - ^ - ^ j • (18) 

These equations imply that once the spacing cjfc between Sk+i and Sk has been decided upon, 
then one can determine £k by inverting the equation at = Aki^k), and then the other parameters 
describing the optimally matched neck can be computed from e^. 

Finally, observe that the matching between ipk{Sk+i) and ipk{.Sk) and the neck defined by the 
choice of parameters above is most optimal in the region of x where the error quantity + 
C(e|/||a;|p) is smallest. It is easy to check that this occurs when ||x|| = 0{e^J^). Hence the radii 
Tk in Step 1 should be chosen as := re^J^ , where the factor of r takes the scaling into account. 

Step 3. Choose dk, and Sk as above. Also, introduce displacement parameters 61,62, - ■ ■ that 
serve to slightly displace the necks from their optimal locations. Define a smooth, monotone cut- 
off function % : [0, cxo) [0,1] which equals one in [0,^] and vanishes outside [0,1]. Define the 
functions : Bi{0) ^ R by 

:= Xdkll/ef )F±,,(6fe, 4 + Sk;x) + (l - xdl^ll • (19) 
Now define the catenoidal interpolation between ipkiSk+i) and tpk{Sk) as 

Nk:={{F+{x),x):£k<\\x\\<£l^^}U{{F-{x),x):ek<\\x\\<el^^}. 
Finally, one can define the finite-length approximate solutions as follows. 

Definition 9. Let K be given. The finite-length surface with parameters a := {ai, . . . , ax} and 
6 := {61, . . . , 5k} is the surface given by 



t^{a,6) := \So\[Bn{p+)UBn{Po)] 
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u 



K-1 



k=l 



u 



U vr 



U Sk\[Bn{p+)uBn{p,)] 
U Sk\[BR{p+)uBR{p,)] 



U 



where vr is the t 



.k=l 
-t reflection. 



Sk \ Br{pj^ 
Sk \ Br{p]^) 



K-l 



k=0 



u 



U ^k'iN, 

U ^k\N, 



U 



fc=0 



3.2 The One-Ended Surface 

The one-ended family of approximate solutions is constructed by attaching a half Delaunay surface 
to a finite-length surface as constructed above. This Delaunay surface has very small necksize 
parameter, so we will need the detailed analysis of these from |9]. 

Step 1. Let K he a large integer. Exactly as in the previous section, construct K perturbed 
spheres of the form Sk for k = 1,...,K and one terminal perturbed sphere of the form S'q 
(which in this case is a normal graph over Sq \ -Bro(Po^))) ^-loiig with perturbed necks V'fe'^(-^fc) 
for k = 0, . . . , K — 1. Then glue these building blocks together using cut-off functions and matched 
asymptotics again as before. This construction should be equipped with the appropriate separation 
and displacement parameters ctq, . . . ax-i and . . . ,5k-i- 



Step 2. The standard Delaunay surface of mean curvature 2 and with the appropriate small neck 
radius can be rescaled by a factor of r and translated along the geodesic 7 until there is overlap 
with the last perturbed sphere Sk of the construction of Step 1. Optimal overlap can be achieved 
because the neck region of this Delaunay surface is to first approximation a standard catenoid. 

To be a bit more precise with this idea, one proceeds as follows. First recall that Delaunay 
surfaces of mean curvature 2 are the surfaces of revolution generated by the functions : M ^ M 
studied in [9l §3] . These functions are periodic with period T and have a local minima at the integer 
multiples of T. Introduce the parameters (which will be fixed once and for all below), and 5k 
and gk (which remain free). Define T := 2 + uk/t and ek '■= Pt{0) and Tk := Ik + r{dK + 5k) 
where Ik is the arc-length parameter for the center of the perturbed sphere Sk- Now parametrize 
a family of translated Delaunay surfaces of mean curvature ^ and period 2r -|- aK via 

H : (t,0) ^ cos(0) , rpT^^^) sin(e) , . 

Note that the t-parameter now corresponds exactly to the arc-length parameter of 7. Define 
truncations of these Delaunay surfaces via 

Dt{aK,5K) := [^{t,e) : t G [Tk-lo, 00) and 9 £ S^} 

where u; is a small number of size 0{re^^'^). 
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Next, one must find the optimally matching Delaunay surface. First note that the neck size ek 
of the Delaunay surface D^{aK, Sr) depends on ax- Now as in Step 2 of the construction of the 
finite-length surface, the equality of the constant and logarithmic terms of the expansion of the first 
neck of D^{aK-, 0) and of the asymptotic region of Sk near determines dx and uniquely in 
terms of the remaining free parameter (tk ■ Now modify (ax , 0) by means of cut-off function 
as in Step 3 of of the construction of the finite-length surface so that it coincides exactly with Sk 
in the region Sk H {{t,x) : ^ruj < \\x\\ < ruj}. Construct also surfaces where 6k 7^ 0, introducing 
mis-match. Denote the family of modified surfaces by D^{aK, Sk)- The definition of the one-ended 
approximate solutions is finally at hand. 

Definition 10. Let K be given. The one-ended surface with parameters a := {gq, . . . ,(Tk} and 
5 := {So,...,Sk} is 



K 



U \ Br{p{) 



k=0 



u 



K-1 



U ^^"'(^^ 



k=0 



UD+{(jk,5k). 



4 Function Spaces and Norms 



The constant mean curvature equation for normal pertTirbations of the approximate solutions 
5) constructed in the previous section will be solved for functions in weighted Holder spaces. 
The weighting will account for the fact that the geometry of S,, (cr, 5) is nearly singular in the small 
neck-size limit that will be considered here. In fact, two slightly different function spaces will be 
introduced. The space C^'°(S^(cj, 6)) will consist of all C^q'° functions on 5) where the rate of 

growth in the neck regions of S^((T, 5) is controlled by the parameter v. The space C^'^(E^^((T, 5)) 
will consist of all C^^'" functions on S^^((t, 5) where the rate of growth in the neck regions of 
S^^((t) is controlled by the parameter u and the rate of growth in the asymptotic region along 
the axis of E^^(cr, 5) is controlled by the parameter v. This latter degree of control is necessitated 
by the non-compactness of S^^(cr, 5). The two types of spaces described above will collectively be 
denoted C*'"(Sr(C) <^)) for brevity whenever needed. 

4.1 Function Spaces and Norms for the Finite-Length Surface 

To begin, one must introduce a number of objects. Define a weight function C,^ : flr{<^,S) — to 
achieve control of the growth of functions on Sr-(C) in the neck regions: 



■t 



r||x|| i^kip) = it,x) G Nk with x G -B^'/2(0) for some k 

Crip) '■= \ Interpolation tpk{p) = it,x) G tpk{Sk) with x G 3^/(0) \ -B^'/2(0) for some k 
r elsewhere 



where the interpolation is such that C,r is smooth and monotone in the region of interpolation, has 
appropriately bounded derivatives, and is invariant under all the symmetries of Sr(<T, (5). For the 
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norms themselves, first introduce the following terminology. If W is any open subset of T>r{a, S) and 
T is any tensor field on U, define 

mo,^:=sup ||T(x)|| and [T^u := sn^ ''"^^"l^fa '^0^^^^^^'' ' 

where the norms and the distance function that appear are taken with respect to the induced metric 
of (f) ) while 3;/ is the parallel transport operator from x to x' with respect to this metric. 
Then, for any function / : W — ^ M define 

k 

\f\k,a,.,U := E ICr^VVlo.W + [C^"-'^VV]a,W . 
i=0 

The norms for the finite-length and one-ended surfaces can now be defined. 

Define a collection of overlapping open subsets of S^(cj, 5) as follows. Let R be a fixed radius 
such that U^, B2ji{p\.) contains all the neck regions of T,^{a, 6) and A := S^(cj, 6) \ [|J^ ^_R(Pfc)] is 
the disjoint union of all the spherical regions of Sr(<7, S). Let Ar := S^((T, d)ri [ (J^ i?2ii(Pfe)\-BR(p^)] 
for any choice of i? € (0, R). 

Definition 11. Let W C {a, S) andueR and a G (0, 1). The norm of a function defined 
on U is given by 

\f\(.k,a.^. ■.= \f\k,a,UnA+ sup \f\k,a,UnAR- 

" -RG(0,R) 

The function spaces that will be used in the case of the finite-length surface arc simply the usual 
spaces C^'"(Sf {a, 6)), but endowed with the C^'" norm. This space will be denoted Ci/'"(S^((T, 6)). 

4.2 Function Spaces and Norms for the One-Ended Surface 

The definition of the weighted norm that will be used in the case of the one-ended approximate 
solution builds upon the norm just defined above. Extend the collection of overlapping open subsets 
used above by defining the points p\. for k > K as the points of 7 upon which the neck regions 
of {uk-i^k) are centered, and then re-defining A and Ar as infinite unions over all fc G N. 
Furthermore, re-define the weight function C,r by leaving it unchanged on E^^((T, 5) \ Df{aKi5K) 
and defining 



r\\x\\ ipkip) = {t,x) G D+{(Tk,5k) with x G 5^y2(0) 

Interpolation ipkip) = G D^^aKjSx) with x G Br/(0) \ Bjii/2{0) 
r elsewhere 



Finally, for any subset U C S^^(fj, 5) and function / G C^^'^(U) define the norm \f\k,a,uM as above. 
A second collection of overlapping open subsets of D^{(Jk,^k) will be now be introduced. Let 
T := be the arc-length parameter corresponding to the point pjv/2 and for any choice oiT >T 
define Ct := S^^(c7, 5) n {(x, i) G x R : i G [T, T + r]}. Also define C := t^^{a, 6) \ Cf. 
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Definition 12. Let hi C S^^((T, 5) and z/, G M and a G (0,1). The C^'p norm of a function 
defined on U is given by 

l/U^'^fl^'l •= \f\k,a,vMr^AnC + sup T""^ l/U.a.i/.Wn^nCr + sup \f\k,a,uMf^ARnCT 

"■■"^ T>f \ iie(o,i?o) / 

The function spaces that will be used in the case of the one-ended surface are the spaces 
Cl'^{to^{a,5)) := {/ G Cf'"(E0-^(a,5)) : \ fLk,c < oo} endowed with the C^'p norm. 

Remark: At first glance, the norm of Definition [T2] looks different from the norm used to study the 
linearized mean curvature operator of near-degenerate Delaunay surfaces in [9, §4]. This is because 
the norm in [HI §4] is defined using a different parametrization (the s-parameter) . However, it is 
straightforward to check, using the estimates of [9l §4] relating the s-parameter to the arc-length 
parameter, that the norm in [9', §4] is equivalent to the Cq'^* norm defined above. 

5 Estimates of the Mean Curvature of the Approximate Solutions 

This section estimates the amount which the approximate solutions T,r{cr,5) deviate from being 
CMC surfaces. This is accomplished by estimating H [Sr (o", 5)] — ^ in the C^'" norm when i/ G (1,2) 
for the finite-length surface and (v, v) G (1, 2) x (—1,0) for the one-ended surface. In the estimates 
that follow, set e := max{e^,efc} and 5 := max{(5fc} and recall that = 0{re^f^). 

5.1 The Estimate for the Finite-Length Surface 

The result for the finite-length surface T,^{a,5) is as follows. 

Proposition 13. Suppose v G (1,2). The mean curvature ofT,^{a,6) satisfies the estimate 



r 



for some constant C independent of r, e, 5 and K. 

Proof. There are several steps: the first two steps are to derive pointwise estimates for the mean 
curvature and second fundamental form of S^((j, 5) with respect to the Euclidean background metric 
in the spherical and neck regions, respectively; the third step is to convert these into pointwise 
estimates for the mean curvature with respect to the actual background metric; and the fourth 
step is to compute the desired C^!_2 norm of F[S^(c7,(5)] - Only the weighted norm will be 
estimated explicitly below since the calculations for the weighted Holder coefficient are very similar. 
Finally, all estimates computed below are independent of K. 

Step 1. The first step it to find a pointwise estimate for B := B [S^(cr, 5)\ and H := If [S^((T, 5)] , 
the Euclidean second fundamental form and mean curvature, respectively, in the spherical region Sk 
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of (cr, 5) . The key to this estimate is to use the formula ([9]) and (jlOp for B and H in terms of the 
graphing function rG and the second fundamental form and mean curvature of in conjunction 
with the estimates from Section [2. 2[ The result is 



l^ll < ^ (l + \G\ + r\\VG\\ + r^llV^Gj 



r 



< -(|G|+r||VG||+r2||v2G||^' 



(20) 



r 

The reason the estimate for |ff — || is so much better is because jC,{G) = in the region being 
considered. 

To proceed with the estimate for |//— it is thus necessary to estimate |G|+r||VG||+r^||V^G||. 
First, in the part of Sk near the points p^^ (here * is A: or /c — 1 as appropriate), where the distance 
to these points can be as small as 0{re^^'^), one can use the Taylor series expansion of G to derive 
the estimate 

\G - ec\og{p)\ + p\\VG\\ + p'^\\V^G\\ < Ce, 

where p := dist(p5^ , •) is the distance function to with respect to the standard, unit-radius, 
induced metric of the sphere. Next, in the part of S/. that is an e-independent distance away from 
these points, G satisfies the estimate 

|G| + ||VG|| + IIV^GII < Ge. 

In both of these estimates, c and G are numerical constants. Therefore near the points one has 



Cre 

\B\\ < — TT- and 
P 



while in the remainder of Sk one has 



G 

\B\\ < — and 
r 



r 



r 



3^2 



< 



Cr^e 



< 



Ge 



2 



Step 2. The next step is to look inside one of the scaled normal coordinate neighbourhoods 
used in the definition of the necks. Again, the calculations are performed with respect to the scaled 
Euclidean metric, and computed in the transition region := Nk^{{t, x) : x £ B 3/4{0)\B 3/4 /„(0)} 
as well as the neck region := n {{t,x) : x £ B 3/4,„(0)} itself. Of course, the estimates in 
Nk are extremely straightforward since N^. is exactly the standard catenoid for which H = and 
||-B|| = \/2£ ll^ll"^- Hence the challenge lies in estimating H and i? in T^,. 

The transition region is the graph of the function F := ijF + (1 — 77) G as in (|19|) . where rj is 
a cut-off whose derivative is supported in ^£^J^ < < e^/^, while F{x) := arccosh(||j;||/efc) -|- 
dk + and G has an asymptotic expansion that matches the asymptotic expansion of F except 
for the mismatch introduced by 5k- Now, the scaled Euclidean second fundamental form and mean 
curvature of a graph are 
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o o ~ 1/2 ~ 

where -D := (l + ||V-F|p) . Write F := u + F where u := {1 — r]){G — F) and observe that 

u = ek6k{l - 7]) + {I - r]){G - F) 

for constants c and C of size 0{e) and functions G and F of size 0(||x|p) + ©(e^HxH"^). When 
ll^ll = 0(e^/^) one has 

u - EkSkil - v) = 0{e^/^) 

whereas ||VF|| = 0{e'^/'^) and HV^FH = 0{e~'^/^). Thus plugging F into the expression for B and 
estimating yields 

II-BII < Ce"^/^ 
Next, plugging F into the expression for H yields 



DH = Au + 



1 



7?.F ./?.. 



'1 + ||VF||2 -^1 + ||VF + Vm||2^ 
1 + IIVF + VuiP 



where AF — (1 + ||VF|p) ^F^iFjF^ij = has been used, which holds since F is the graphing 
function for the catenoid which has zero mean curvature. Therefore one can estimate 



\H - 2 + ekSkAr]\ <C 



where C is independent of e. 



Step 3. The next step is to compute the pointwise norm of H := H[T,^ {a,6)] with respect to 
the actual background metric. Lemma [2] relates H to H and B, and the form of the relationship 
that is germane to the current derivation is expressed most succinctly in equation (jl2p . Thus to 
proceed, one must substitute the estimates from Steps 1 and 2 into this formula and estimate the 
curvature terms as in Section 12. 3i Again, this should be done in the different regions identified 
above. Consider first the spherical region Sk, where ||y|| = 0{r) and 



H 



< 



< 



H 



H 



+ |i?(0)|||i?|| + |i?(0)|| 

+ C7f||y|| + ||y||2||s||) 
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< < 



c 



C - + r 
V ^ 

near and away from p^, respectively. 

Consider now the neck and transition regions Ni^. First, by reversing the scahng used in Step 2, 
the second fundamental form and mean curvature in N]^^ measured with respect to the Euclidean 
metric, satisfy the estimates 



a \-ek6kAr] 



< 



C 



and 



I-BII < 



C 



re 



1/2 • 



Now these estimates can be plugged into the expansion of the mean curvature as above, except that 
Y is now the position vector field of ip~^{Nk) relative to the center of the normal coordinate chart 
used in the construction of the neck — namely the point p^. Hence ||y|| is uniformly bounded by 
max{r|G|, r|F|} < 0{re\ log(e)|). With this in mind, one obtains the estimate 



H 



< < 



C 



-(1 + r^e\ log(e)| + r'^e^/^\ log(e)|2) 



1 + rh\ log(e)| + ""'^'1/°!,^^^'' + reA\A7j\ ] . 



\ r \ \\x\ 
in the transition region of the neck and in the neck itself, respectively. 
Remark: In all of these estimates, the TC term from Lemma [2] is always negligible. 



Step 4. The remaining task is to estimate \H — ^\ in the norm for u G (1, 2). This estimate 
will be derived by estimating the supremum of Cr"'' I -f^ - - 1 in the three regions of S^(cr, 6) identified 
in the previous steps. For now, the only assumption that will be made about e is that e <C r. First, 
consider a spherical region Sk away from the points p^. In this region (^r = r so that 

2 



H 



< Cr 



3-u 



Also, in a spherical region Sk near the points p^,, one has instead Cr ~ dist(p' 

2 



rp. Thus 



H 



< sup Cirpy-" r + 

pG[£3/4,R'] V P' 

Next, in the neck region one has Crix) = r\\x\\ in the local coordinates used to define N^- Hence 



H 



< sup Cr 

£<||z||<e3/4 



1-.,|^||2~. 



1 + 



+ r'^e\log{e)\ + 



r^e^l log(e) 



+ ekSk\Ar]{x) 



< Cri-^((1 + r^e\ log(e)|)e3(2-)/4 + ^ ^2^3-.| ^^^^^^^2 ^ ^^^-1/2^ 
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Therefore consolidating these three estimates yields 



H 



Cl-2 



< C max |r 



■i-u ^l-,.^3/2-3i//4 ^_^l-i.^l-3;./4 



e 



which is the desired weighted supremum norm estimate. The full estimate follows once the 

estimate for the Holder coefficient has been computed. As indicated above, this computation is 
more involved but very similar, and yields the same result. □ 



5.2 The Estimate for the One-Ended Surface 

The estimate of the mean curvature of the one-ended surface $]^^(cr, (5) is a very straightforward 
extension of the results obtained in the previous section. In fact, the result is the same. 

Proposition 14. Suppose v G (1,2) and v G (—1,0) is sufficiently close to zero. The mean 
curvature o/S^^(ct, 5) satisfies the estimate 
2 



H 



< Cmaxjr^--, ^5-.^l/2-3./4^ ^l-.^3/2-3./4^ ^^1-.^1^3./4| 



for some constant C independent of r, e, 5 and K. 

Proof. The computations of Proposition [13] give the estimate of the mean curvature of the part 
5]^^(cT, 5) constructed from spheres and catenoids. It thus remains only to compute the estimate 
of the mean curvature of D^{aK,SK) (actually, the part of D^^gki^k) that is an un-perturbed 
Delaunay surface). Once again, the key is to use Lemma [T] in the form of equation (|12|) . but this 
time realizing that the mean curvature of [gk^^k) with respect to the Euclidean background 
metric in a tubular neighbourhood of 7 is exactly equal to ^. Let p be any point on 7 and let p' be 
any point on D^{aK-, Sk) H Br{p). Then by equation (fT2]) . 

2 



H{p') 



<C\R{p)\{\\Y{p,p')f\\B{p')\\ + \\Y{p,p') 



where R is an expression that is linear in the components of the ambient Riemannian curvature at 
p and Y{p,p') is the position vector field of with respect to p. The constant C is independent of r 
and e. Since the curvature is exponentially decaying along 7, the same estimates from Proposition 
[T3] continue to hold and yield the estimate needed here. □ 



6 The Solution up to Finite-Dimensional Error 

6.1 The Finite-Length Surface 
6.1.1 Strategy 

Let fi : C^'"(S^((T, 5)) Emb{Ti^ {a, 6), M) be the exponential map of S,f(cj, 5) in the direc- 
tion of the unit normal vector field of S^((T, 5) with respect to the backgroung metric g. Hence 
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firfi^ri^i^)) is the scaled normal deformation of T,r{(^,S) generated by / S C^'"(S^(cr, (5)). The 
equation 

H[^^,f{t^{a,6))]=l (21) 

selects / G Cv'"(5i^(cr, (5)) so that jj,rf{Ti^{a,6)) has constant mean curvature equal to |. In 
addition, the function / will be assumed symmetrical with respect to all the symmetries satisfied 
by S^((T, 5). Using a fixed-point argument together with a suitable choice of weight parameters, 
the equation (j2ip will be solved up to a finite- dimensional error term. This means that a solution 
of 

H[firf{t^{a,6))] = l+£ 

will be found, where £ belongs to a finite-dimensional subspace of functions that will be denoted 
and specified below. At first glance, the error £ will come from terms in the solution procedure 
that are not sufficiently small. However, the true reason for the presence of £ is geometric and will 
be explained in Section [TJ where we show how to eliminate it. 

To begin, write ff[/Ur/(S^((T, 5))] = HlTiric, 6)]+C{r f) + Q{r f)+7i{r f) where £ is the linearized 
mean curvature operator, Q is the quadratic remainder part of the mean curvature and TC is the 
small error term as in equation (jlip . The first step is to construct a suitably bounded parametrix 
TZ : C'l^^gl^r'l^i <^)) ~^ Cu'^{Ti!^ {a, 6)) satisfying C olZ = Id + £ where £ maps into VV^. Now the 
Ansatz f := ^TZ(^w — //[S^(cj, 5)] + ^) transforms the equation (|21|) into the fixed-point problem 

w = -QoTz(^w - H[t^{a,6)] + - H oTZ (^w - H[t^ {a,6)] + . (22) 

up to the finite-dimensional error term in . The remaining task is to show that the mapping 
A^r : C'^'^2(^r {(^^ ^)) ~^ C^-2{^r ^)) given by the right hand side of ([2^ is a contraction mapping 
onto a neighbourhood of zero containing [S^(cr, (5)] — ^. Once this is done, then one has solved 
the equation ([26|) up to a term in . 

6.1.2 The Linear Analysis 

We now find a parametrix TZ satisfying C olZ = id + £ where £ has finite rank. In each normal 
coordinate chart, C is close to the linearized mean curvature operator with respect to the Euclidean 
metric C{f) := Af + Patching together inverses for the latter operator gives a parametrix 

with an error which decomposes into terms which are genuinely small, and those which together 
constitute £. 

First, denote Ann-rip) •= ^rip) \ Bt/2{p) ^-iid define subsets of T,^{a,6) by 

■.= t^{a,5)\[Br{pi)UBr{pi_^)] 
Tf;;^:=t^{a,6)nAnnr{pi) 
Tk~ ■= (cj, 6) n Annr{pl_i) 
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Now define the smooth, monotone cut-off functions 



Xneck,k \ 



Interpolation 


1 

Interpolation 




elsewhere 



xeS^ 
elsewhere 



so that Y^f^ x^ext A: + X^fc Xneck k ~ ^ cut-off functions are invariant with respect to 

all symmetries satisfied by (a, 6) . 



Proposition 15. Let u E (1,2). There is an operator 7^ : C^'^2i^r (f^><^)) ^ Cu"{T,^ {a,S)) that 



satisfies Con = id-£ where S : C°!?2(Sf {a, 5)) ^ . Here is a finite- dimensional space 
that will be defined below. The estimates satisfied by TZ and S are 

for all w G C°:^2(^f (^^^ where C is a constant independent of r, e and S. 

Proof. Let w G C^'^2{'^r ^)) be given. The task at hand is to solve the equation C{u) = w+£{w) 
for a function u G C^'"(I]^((T, (5)) and an error term £{w) G VV^. To begin, introduce four radii 
Ti < T2 < Ts < T4 <^ r with the property that the supports of the gradients of the cut-off functions 
X? and x* do not overlap for i ^ j. These radii will need to be further specified; and this will be 
done in the course of the proof below. 

Step 1. Let Wneck,k '■= "^X^eckk- '^^^^ function has compact support in J\f^ and thus can be 
viewed as a function of compact support on the standard scaled catenoid. Now consider the equation 
JO-n{u) = Wneck,k the Standard scaled catenoid, where £jv is the linearized mean curvature 

operator of the scaled catenoid re^N with respect to the Euclidean metric. Then the pull-back 
of C to the standard scaled catenoid is a small perturbation of Cn. By the theory of the Laplace 
operator on asymptotically flat manifolds, the operator /Z^v is surjective onto C|^^2(^^fc-^) when 
v G (1,2). Hence there is a solution Uneck,k as desired, satisfying the estimate \uneck,k\c'^-'^ — 



C\w\ 



where these norms can be taken as the pull-backs of the weighted Holder norms being 



used to measure functions on (cr, (5). Finally, extend Uneck,k to all of E^(a", 5) by the definition 

'^neck,k ■ ^ neck, kX neck, k "^^^ Set Uneck • Y2k^'"-^ck,k- 
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Step 2. Define Wext,k '■= {w — ^{uneck))x^ext k- Then Wext,k is a function of compact support on 
and thus can be viewed as a function of compact support on the sphere Sk \ {p^,P^}- Now 
consider the equation Cs{u) = Wext,k, where Cs is the hnearized mean curvature operator of the 
sphere of radius r with respect to the Euchdean metric. Then the pull-back of C to the sphere is 
a small perturbation of Cs- It is not a priori possible to solve the equation Cs{uext,k) = Wext,k on 
Sk because of the one-dimensional kernel of Cs that is invariant with respect to the symmetries 
that have been imposed. A basis for the kernel is given by the function J/t : 5^ ^ R defined 
by Jk ■= di^j §i) where g is the Euclidean metric and N is the unit normal vector of Sk- Let 
'^^xt k •~ '^ext,k — {wext,k, Jk)Jk be the projection ofwext,k to the orthogonal complement of spanjjfc} 
with respect to the Euclidean L^-inner product denoted by (•,•). Now there is a solution of the 
equation Cs{uext,k) = w^xt,k satisfying the estimate \uext,k\c2,cc < C\ 

'^'ext,k\c'^''^ where | • \Qk,a is the 

scale-invariant C^'" norm where derivatives and the Holder coefficient are weighted by appropriate 



u-2\ 



W 



factors of r. Furthermore, one has |u^ext fclco-" — ^l^ext,fc|^o,a < 

A modified solution satisfying weighted estimates can be obtained as follows. First, write 
Uext,k ■= Vext,k + o-t^k + H% where Vext,k satisfies \vext,k{^)\ < CT2~'^dist{x,pl)'^\w\^o,^^ near p^, 
while the function r]^ equals one near and vanishes a small but e-independent distance away 
from these points, and G R satisfies \a^\ < CT2~'^\w\^o,a . This decomposition is achieved by 
studying the Taylor series expansion of Ue^i^fc near p^ and using the symmetries satisfied by Uext,k- 
Finally, by adding the correct multiple of Jk to Uext,k on each 5^, one can arrange to have = 
for every k > 1 and = by symmetry. To extend the functions Ugxt.k to all of S^(cr, 6), define 
Uext ■= Vext + A where Vext ■= Yjk x2:t,k'^ext,k and A := ".t^k'X.neck.k- Then one has the estimate 

|f}ea:tlri2.a + To '^\A\„2.a < Cllt^l^O.a . 

step 3. Let u(^) := Vext + Umck and £^^\w) := -Y^kX2xt,ki^ext,k,Jk)Jk - Y.kX2xt,k^s{A). By 
collecting the estimates from Steps 1 and 2, one has I The claim is that 

\C{u^^^) -w - £'^'^\vj)\^o,c^<9\w\^o,o.^ and \£^^\vj)\^o,c < Ct^~'^\vj\^o,c^ (23) 

where 9 can be made as small as desired by adjusting ri, . . . ,r4 and e suitably. The consequence 
is that one can iterate Steps 1 and 2 to construct sequences n*^"^ and £^^\w) that converge to 
u := TZ{w) and £{w) respectively, satisfying the desired bounds. 

Therefore to complete the proof of the proposition, it remains to compute the estimates given 
inl23[ The idea is to exploit the fact that C differs very little from Cn and >Cs in the regions where 
u^^^ equals Vext and Umck-, all while taking into account the effects of the cut-off functions. With 
this in mind, and using the notation [£, x\{u) := C{xu) — xC{u), one has 

£(n(^)) = (£ - Cs)iVext) + '^[Cs,xll:t,k'\(^e^t,k) + '^Xllt,k^s{Vext,k) + C{Uneck) 

k k 

= {C- Cs){Vext) + ^[Cs,X"'ext,k\{Vext,k) + J2 Xllt,k'^ ext,k 
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- Yl XZt,ki^ext,k, Jk)Jk - Xllt,k^s{A) + C{Uneck) 

k k 

= {C- Cs){Vext) + X]['^-S''XlLt,fc](^e^t,fc) + Z]^ext,fc(^ - ^Uneck)) 
k k 

-Y^^ext,k{'^ext,k,Jk)Jk - Y^''ext,k^s{A) + C{Uneck) 
k k 

= {C- Cs){Vext) + Xll-^S, XlLt,J(^ext,fc) + Yl + Y Xneck,k^(^neck) 

k k k 

k k 

= {C- Cs){Vext) + J^[/:S,XlLi,fc](^exi,fc) +Y^neck,ki^ " ^N){Uneck) 
k k 

- Y ^iit,k{'^ext,k, Jk)Jk - Y + w 

k k 
since C{Uneck) = - C,N){Uneck) + Ylkl^N ,x''neck,k\^'^neck,k) + T.k X^neck,k^N iUneck,k) ■ The factS 

xllt,kX7^t,k = xZt,k and Xneck,kXneck,kXneck,k = xZck,k ^nd Xneck^^N , xZckJ = ^ ^ave also been 
used in the calculations above. The result now follows since both C — Cs and C — Cn can be 
handled using the estimate p6p while 

I [Cs, Xllt^k]^extAcl'-, ^ ^l^-t.'^lc^'-CsuppCVxS,,,)) ^ ^ (^) \'"\cl'-^ 



can be made as small as desired by adjusting the ratio ti/t2. 

The iteration leading to the exact solution of the equation C{u) = w + £{w) now works as 
follows. Using the steps above, for every n > one has functions n*^"^ G c3'°(S^((T, 5)) and 
^(„) ^ c^'^^it^ {a, 5)) satisfying = u'("-^) + £:(n("-i)) + u-^") and w^o) = w along with the 

estimates 

l-u^^^La,. + |£:(u;("-i))L2,. < C|u;("-^)Lo,. and ju-^^^Lo,. < -\w'^'^-^^\^o,c . 

Consequently the series u := — l)"'''^'"^") and £{w) := l)"'"''^'? (^f ^"~^^) converge in 

the appropriate norms and satisfy C{u) = w + £{w) along with the desired estimates. □ 

The definition of the finite-dimensional image of the map £ : C^'^2{^r i^-:^)) ~^ is a by- 
product of Step 3 of the previous proof. 

Definition 16. Define 

:= span{x:^,,, J. , x7xt,kCsi4x:L,,,) : k = 1, . . . , K - 1} U {xHt^K^K} • 

6.1.3 The Non-Linear Estimates 

The next task is to find estimates for the C^'" norm of the quadratic remainder term Q and the 
error term 7i for the finite-length surface. We do this by combining the estimates from Section [2] 
with the specifics of the construction of S^((T, 6) from Section [3l 
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Lemma 17. Pick x G S^((t, 5). Then x belongs to one of the normal coordinate charts used in the 
construction of T,^ (u, 5) where the second fundamental form with respect to the Euclidean metric 
is B{x). At this point, the estimate 

r{\\B{x)\\\f{x)\ + \\Vf{x)\\)<Cr^\f\^.,. 

holds for all f G C^'"(S,^(cr, (5)), where C is a constant independent of r, e and 5. 

Proof. Collecting the estimates for 13 from Section [5l one finds that 

Cr{x)\\B{x)\\+Ch^)WB{x)\\<C (24) 

where C is some constant independent of r, no matter where x is located in T,^{cr, 6). Consequently, 

\\B{x)\\\f{x)\ + ||V/(x)|| = Cr\x){\\Bix)\\Cr{x) • C-"(x)|/(x)| + C^'^(x)||V/(x)||) 

which yields the desired estimate. □ 

Hence it is true that r(|/|||i?|| + ||V/||) <^ 1 can be ensured by keeping |/|^2,a small enough. This 
condition validates the computations of Section [2l The following estimates are a consequence. 

Proposition 18. There exists M > so that /2 G C^'"(S^((j, 6)) for v G (1, 2) and satisfying 
\fi\u2.a + 1/2 1 (-,2, a < M, then the quadratic remainder term Q satisfies the estimate 

|Q(/i) - Q{f2)\cO.^^ < Cr''~'^\fi - /2|(^2,« max{|/i|(^2,«, |/2|^2,a} 

where C is a constant independent of r, e and 6. 

Proof. Choose a point x G T,^{a,6). Then x belongs to one of the normal coordinate charts used 
in the construction of S^(cr, 5). If r(|/j|||i?|| + ||V/i||) is sufficiently small, it is possible to invoke 
Lemma [7] along with the estimate of the previous lemma and immediately deduce 

CMQ(/i) - Q(/2)| < CCrl/i - /2I • max {\fi\\\Bf + \\Vfi\\\\Bf + \\V fi\\\\VB\\ 

+ iivv.iiPii + i/d + (1 + ni)iiv/.ii) • C? 

+ CCM|V/i - V/2II • max {\fi\\\Bf + \\Vfi\\\\B\\ + \fi\\\VB\\ 

+ llvVill + (1 + \\y\\)\fi\ + (1 + rf )I|V/.|| + ||V/i||||VVi||) • Cr 
+ CCMlVVi - VV2II • max {\fi\\\B\\ + \\Vfif) 
< c(r'^-i +r2-2 +r'^+i +r'^+2)|/, _ /2|^2,. max {|/,|^2..} 

since 1/C < HyHC^"*^ < C for some universal constant C. The desired estimate follows. □ 
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Proposition 19. There exists M > so that if fi, f2 S Cu'°' {T,^ {a, 6)) for v G (1, 2) and satisfying 
|/i|^2,a + 1/2 1 ^2, a < M, then the error term TC satisfies the estimate 

\n{fi) - n{f2)\co,^ < CrVi - Mc^-" 

where C is a constant independent of r, e and 6. 

Proof. Similar computations. □ 
6.1.4 The Fixed-Point Argument 

We are now in a position to solve the CMC equation up to a finite-dimensional error. Let 

E := H[t^{a,6)] and R{r,e, 6) := max {r^'" , r^-^ e^l'^-^^l^ ^ ^i~u^3/2-3u/4^^^^i~3u/4y 

tionally, assume r^ < e < r^ 1 and 6 < e^l"^ . All of this is justified a posteriori. The following 
estimates have been established. 

• The mean curvature satisfies li^l^o,!! < CR(r,e,5). 

• There is a parametrix IZ satisfying C oTZ = id — £ where £ maps into the finite-dimensional 
space VV and \7l{w)\^2,a + \£{w)\^2,a < C\w\fjO,a for all w E C|^^2(^r(<7, (5)). 

• The quadratic remainder satisfies |Q(/i) — Q{f2)\(jO,a < Cr'^^^\fi — /2|^2,q maxj { |/j|p2,a } 
for all /i,/2 G C^'"(I]^(cT, 5)) with sufficiently small C^'" norm. 

• The error term satisfies \TC{fi) — 'H{f2)\ , < CrVi - /2lc2.- for all /i, /2 G ^''"(Sf (a, 6)) 
with sufficiently small C^'" norm. 

One can now assert the following. 

Proposition 20. There exists w := Wr{cr,5) € C'^'"2{^r i^^ ^)) ^'^^ corresponding f := fr{cr,5) € 
{a, 6)) defined by f := ^n{w - H[t^ {a, 6)] + so that 

H[f,rf{t^{a,d))]-^ = -£ (25) 

where £ € . The estimate \ f\^2,a < Cr^^R{r,£,6) holds for the function f , where the constant 
C is independent of r, e and 5. Finally, the mapping (cr, (^) ^ f{a,5) is smooth in the sense of 
Banach spaces. 

Proof. By the last three bullet points above, the map w 1— > Mr{w) satisfies 

\Nriwi) - Nr{w2)\r^Q,c. < C (r""'^ R(r, E, 5) +r'^)\wi - W2Uo,<. 

^i/-2 ' ^1^-2 

where C is independent of r and e. Since r^^^Rir, e, 5) and can be made as small as desired by 
a sufficiently small choice of r and e with < e < and 5 < e^^^, it is thus true that Mr is a 
contraction mapping on the ball of radius R{r,£, 6) for such r, e and 6. Hence a solution of (|25|) 
satisfying the desired estimate can be found. The dependence of this solution on the parameters 
(cj, 6) is smooth as a natural consequence of the fixed-point process. □ 
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6.2 The One-Ended Surface 



6.2.1 Strategy 

The strategy for solving the CMC equation ()2ip in the case of the one-ended surface (cr, 5) 
must be modified in order to take the non-compactness of S^^(cr, (5) into account. In fact, 
the modification required can be understood by considering the outcome of the hnear analysis, 
specifically the nature of the parametrix for C. In this case, the outcome of the construction 
the parametrix, which will mimic Proposition [15] as closely as possible, will be a parametrix 
TZ : C'^'\j,{toi^{a,d)) CJ^(SO^(cj, J)) © V satisfying CoTZ = id + £. The operator £ again 
maps into a finite-dimensional subspace W'^^ . The subspace V is the new ingredient, and can be 
explained as follows. First, let Jfj^i for s = 0, 1 be the bounded and linearly growing Jacobi fields 
of the standard Delaunay surface and define the space 

V := Span{x£)el'^i>ei ' XDel^Del} 

where X^el ^ smooth, monotone cut-off function that transitions from zero to one in the neck 
region where the Delaunay end of Ti^^{a,6) is attached to the finite part of Tj^^{a,6). The 
reason V is needed is simply because C : C^l^{t^^{a,6)) ^ C°'^2^{t^^{a,6)) is not surjective 
but becomes so when growth like the first non-decaying Jacobi fields of C is permitted. But now 
the fact that one component of the solution of the linearized problem does not decay forces the 
modified approach that will be outlined in the next two paragraphs, since the quadratic remainder 
of the mean curvature will not behave appropriately for this component. An approach similar to 
the one proposed below has been used in [17]. 

To compensate for the non-decaying component of the solution of the linearized equation, one 
proceeds as follows. Let TZ^^^ denote the component of TZ mapping into Cl'"{tOE{a,5)) and let 
TZ^"^^ be the component of TZ mapping into V. Furthermore, if TZ^'^\w) = diXhei'^Dei ~^ '^'^Xhei'^Dei^ 
then TZ^'^\w) := {ai{w),a2{w)) despite the slight abuse of notation that this represents. Now the 
equation that needs to be solved is still 

H[^irf{t?^{a,6))] = '^ (26) 

but for / S pi^r^ i^^ ^))- Recall that the last two free parameters of T,^^{a,5), namely aK 
and 6k, parametrize asymptotically non-trivial deformations of S^^(cr, (5). Namely, these cause 
the period and location of the entire Delaunay end to change. The idea for converting ()26p into a 
fixed-point problem that can be solved in the standard way, is to associate 7^(l) with / and Tl^^) 
with the parameters {(7k,Sk) in an appropriate way. 

This can be done as follows. Recall that there are specific values ax = ctk and 5k = which 
produce optimal matching in the assembly of T,^^{a,6). With slight abuse of notation, write 
T,^^{ai,a2) := I]^^(o", J) with {(Tk,5k) = {<^K + 01,02)- Given the flexibility one has in the choice 
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of V, one can arrange that 

^"-^ a=0 

for i = 1,2. Consequently, the Ansdtze 

(ai,a2) := ^TZ'^'^ {w - H[f,rf{t?^{a,5))] + f) 

along with the expansion of the mean curvature into its constant, linear and quadratic and higher 
parts transforms the equation (j26p into the fixed-point problem 



«^ = -Qao7^« (^w-H[firf{t^''{a,6))] -non(^w-H[firf{t^''{c7,6))] + '^^ (27) 

up to a term in VV*^^, where Qa denotes the quadratic remainder of the mean curvature of 
S,^^ (01,02). One must now show that the mapping Mr : pi'^r'^ i'^^ ^)) ~^ ^^^^i^-,^) given 

by the right hand side of (j27p is a contraction mapping onto a neighbourhood of zero containing 
[S^^(cr, (5)] — J,. If so, then one has solved the equation (126p up to a term in VV*^^. This finite- 
dimensional error term must of course still be dealt with in order to find a true CMC surface near 
to S^^((T, 5). This will also be carried out in Section [71 

6.2.2 The Linear Analysis 

To begin the construction of the parametrix in the case of the one-ended surface, one must first 
define an additional set of partitions of unity for S,^^(cj, 5) as follows. Denote Cyljf. := {{x,t) € 
M : \\x\\ < r and t > f} for any r S M and then define := Ti^^{a,6) H Cyl^. as well as a 
smooth, monotone cut-off function xhei that equals one in 2?"^ and vanishes in Sj9^(o", 5) \ 'D'^~^ . 

A second important ingredient that will be used in the construction of the parametrix is a 
careful analysis of the properties of the linearized mean curvature operator of a near-degenerate 
Delaunay surface with respect to the C^'p norm. This was carried out in [9] and the relevant results 
from §4] can be adapted to the needs of this paper and will be quoted whenever they are used 
in the proof given below. 

Proposition 21. Let £ (1,2) x (-1,0). There is an operator TZ : C'^'^2^^^{t,^^ {a,5)) 

Cj"(E^o^(a, 6))®V that satisfies £on = id-£ where £ : C°'\^{tOE{a, 6)) ^W. Here W^E 
a finite- dimensional space that will he defined below. The estimates satisfied by TZ and £ are 

\TZ{w)\^,2,^^^ + \£{w)\^,2,^ < C\w\^o,^^_ 
for all w G C^'^^^C^^^ i(^^^))> where C is a constant independent of r, e, 6 and K. 
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Proof. Let w G C^^2p(^r ('^''^)) be given. The solution of the equation C{u) = w + £{11}) will 
be constructed broadly along the same lines as in Proposition [15] in that local solutions on the 
spherical constituents, the necks and the Delaunay end of T,^^{(t,6) will be patched together. In 
this case, however, a preliminary step is needed to reduce the interaction between the Delaunay 
end and the finite part of T,^^ {a,6). This amounts to showing that one can assume that w = in 
a large part of Tj^^{a, 6). 

Step 0. Define Wmid '■= Xoei ~ Xd^i)w. Using the methods of Proposition [151 one can solve 
the Dirichlet problem C{umid) = Wmid + £{wmid) and u^id = on d[V^''^ \2?*2A'j _ ^j^he estimate one 
obtains is lumidlr'^.'^ + \£(wmid)\r<2,a < Clwmidlr^o.a < Clwl^o.a where C is independent of r, e 

'-'0 ^v^2 '-'v-2,i> 

and K. Then Umid can be extended to all of S^^((T, 5) by defining Umid '■= xEeX^ ~ XrHf^) '^mid- 
If one now solves C{u) = w — C{umid), then the function u + Umid solves £,{u + Umid) = w. The 
advantage will be that the function w — C{umid) vanishes in 2?*^ \ I?*2k ]-,^^ g^jj^ satisfies the same 
estimates as did w. One can also assume that K can be as large as desired. 

Step 1. Let w E Cl'"^^^{t'^^ {a, 5)) be given and assume that w = Qui V^'< \ P*2^. Consider the 
equation C{ufin) = Wfin + £{wfin) where Wfin := (1 — X^Dd) ^ view the function on the right 
hand side as being defined on the finite-length surface. Using the methods of Proposition 1151 one 
can find a solution Ufin G C3'"(S^(ct, 5)) satisfying |n^„|p2,a + <S(w^„)|pO,a < C\w\^o,a for some 
constant C independent of r, e and K. Note that one must view 

ptif \ P*2K as being very close 

to a union of spherical and neck regions in order to do so. Furthermore, by carefully analyzing 
the iteration process that leads to the solution, the fact that supp{wfin) C Ti^^{a,6) \ P*^ implies 
that for X G T,^^{a,6) near (9D*2if one has {x)\ufin{x)\ < Ce~^\w\^o,a for some constant C 

u — 2,1' 

independent of r, e and K. Extend Ufin to all of Tj^^{a, 5) by defining Ufin := (1 — Xnei^^) '^fin- 
Step 2. For w € C°^2 ui^r'^i^^ ^)) fro™ Step 1, set w^ei '■= X-qIi ^ and consider now the equation 
C^oiy^Dei) = 'W^De/ but vicwed as an equation on the complete Delaunay surface D^.- Then using the 
methods of [HI §4] and [15] along with the condition v G (—1, 0), there is a solution of this equation 
which can be decomposed as umi = voei + a^idXDdJhei + "-DdXDdJDd- Here vod S Cl'^{Dr) 
and one has the estimate \v]:iei\^2,a + \a\)f,i\ + \o^Dei\ — C\w£)ei\fjO,a for some constant indepen- 
dent of r and e. Furthermore, because supp{w£)ei) ^ Tj^^{a,d) n P*^ one can arrange to have 
{x)\ufin{x)\ < Ce~^\w\^o,a for x G T,^^{a,6) near dV^^^ for some constant C independent of 
r, e and K. Now it is possible to view Ti^^{a,S) H 2?*^ as a graph over Dr H P*^ and hence u^ei 
can be viewed clS Si function defined on 6) D P**'. Extend this function to all of S^^(<t, 6) 

by defining uod ■= X%i'' UDd- 

Step 3. The estimate of C{u)—w up to a finite-dimensional error term proceeds in the same way as 
in Step 3 of Proposition [151 The extra exponential decay of ujin and UDd near \P*2A'j ensures 
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that the cut-off errors that arise there are small. Consequently one can iterate the steps above and 
find the desired solution TZ^w) G C^'"(I]^^(cr, 6)) © V and satisfying the desired estimate. □ 

The definition of the finite-dimensional image of the map £ : C'^'^2i'^r^i^j^)) ~^ VV^^ is once 
again a by-product of the previous proof. 

Definition 22. Define 

W^^ := span{x-,,4 , x7.t,k^s{vtxZck,k) : k = 0, . . . , K - 1} U {xH^^kJk} ■ 

6.2.3 Tlie Non-Linear Estimates 

The estimates for the u norm of the quadratic remainder term Qa in the case of the one-ended 
surface are very similar to the analogous estimates for the finite-length surface. First, the result of 
Lemma [T71 continues to hold because the calculations are essentially identical, the only difference 
being the need to multiply by factors of e"*^^ along the end of I]^^(fT, S). These growing factors are 
compensated for by the exponential decay assumed for the function /. Consequently it is possible 
to make ||i?|||/| -|- ||V/|| pointwise small everywhere by choosing |/|^2,c. sufficiently small. 

Next, the non-linear estimate analogous to Proposition 1181 follows similarly because the terms in 
Qa and H coming from the background metric decay exponentially. One has the following results. 

Proposition 23. There exists M > so that if i\,h G Cl'^i^'^^ for v £ (1,2) and 
v G ( — 1,0) sufficiently close to zero and satisfying |/i|^2,a -|- |/2|^2,q < M, then the quadratic 
remainder term Qa satisfies the estimate 

\Qa{fi) - Qa{f2)\c0:^ < r^-^CI/i -/2|^2,.max{|/i|^2,.} 

where C is a constant independent of r, e and 6. 

Proposition 24. There exists M > so that if fi,f2 G C^;" (S^(cr, 5)) for v £ (1,2) and v e 
(—1,0) sufficiently close to zero and satisfying |/i|^2,a -|- |/2|^2,c« < M, then the error term 7i 
satisfies the estimate 

|W(/l)-W(/2)|^0,. <CrVl-/2|t;2.. 

where C is a constant independent of r, e and 6. 

6.2.4 The Fixed-Point Argument 

The fixed-point argument in the case of the one-ended surface is again broadly similar to the ar- 
gument in the case of the finite-length surface. However, the strategy adopted for dealing with 
the non-decaying component of the parametrix requires some additional care. As before, let 
E := H[tOE^a,6)] - f and R{r,e,6) := maxjr^"'^, r5-i.£i/2-3v/4^ ^i-u^3/2~3u/A^^^i^u^i-3u/4y 

Additionally, assume r'^ < e < r^ <^ 1 and 6 < e^l'^ as before. The following have been established. 
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• The mean curvature satisfies |£'|(^o,a < CR{r,e,5). 

• Tliere is a parametrix TZ satisfying C olZ= id — 8 where £ maps into the finite-dimensional 
space W and \£{w)\^2,a < C\w\^o for all w E C^'^^ pi^?'^ i'^^^))- But now TZ decomposes 
as 7^(l) +71(2) g^^^ |7^{l)(w;)L2,. + ||7^(2)(u;)||- < C\w\^o,c. for all w £ C'^'^^ pi^?^ i^^^ ^))- 

• The quadratic remainder satisfies I Qa(/i) — Qa(/2)Uo,a < Cr*^"^!/! — /2|^2,q maxj | |/iL2,a | 
for all /i, /2 G C^^'p 5)) with sufficiently small C^'p norm. 

• The error term satisfies |7^(/i) — 7^(/2)Uo,q < Cr^|/i — /2U2,q for all/1,/2 G C^'p (Sf (cr, 5)) 
with sufficiently small C^'p norm. 

One can now assert the following. It's proof is entirely analogous to the proof of Proposition I2UI 

Proposition 25. There exists w := Wr{cr,6) G ui^?^ ^)) > corresponding f := fr{a,6) G 

Cj?(S^0^((T,5)) and (01,02) G V de/^ned &y / := ^TZ^^^w - H[t^^{a,6)] + f) and (01,02) := 
7^(2)(y; _ H[toi^{a,6)] + so i/iai 

ii-[//,/(S,o^(oi,02))] - ^ = -'? 

where £ G W*^^. T/ie estimate \f\^2,a < Cr^^R{r, e, 6) holds for the function f, where the constant 
C is independent of r, e and 5. Finally, the mapping (cr, 5) 1-^ f{a,6) is smooth in the sense of 
Banach spaces. 

7 Solving the Finite-Dimensional Problem 

7.1 Strategy 

It has now been established for both families of surfaces that if r, a and 5 are sufficiently small, 
one can find Wr-((7, 5) G C*'"(Xlr(c, (5)) and corresponding fr{cr,S) G C*'"(Sr.(c7, (5)) so that 

2 

where £r (c, 5) is an error term belonging to the finite-dimensional space W* and whose dependence 
on the free parameters in Sr(o", 5) has been indicated explicitly. To complete the proof of the 
main theorem, we must show that it is possible to find a solution where these error terms vanish 
identically. 

Consider the balancing map defined by 

Br{a,6) ■.= TT{£r{a,5)) (28) 

where vr : VV* — > R"' is a suitable bounded projection operator, where d is the dimension of VV*. 
(The operator n will be a certain bijective L^-orthogonal projection onto a finite-dimensional vector 
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space.) Note that Bj. is a smooth map between finite-dimensional vector spaces by virtue of the fact 
that the dependence of the solution fr{a,5) on (a, 5) is smooth and the mean curvature operator 
is a smooth map of the Banach spaces upon which it is defined. It will be shown using the implicit 
function theorem for finite-dimensional maps that for every sufficiently small r > 0, there exists 
{(7,5) := {a{r), 5{r)) for which Br{cr,5) = identically. It is at this stage that the precise nature 
of the scalar curvature of the background manifold M enters the picture: the behaviour of the 
scalar curvature along the geodesic 7 enters into the selection of the parameters a and 5 to first 
approximation . 

7.2 The Balancing Formula 

The projection operators that will be used to study the finite-dimensional error £{wr{cr, 5)) in the 
case of the finite-length surface and in the case of the one-ended surface can be defined as follows. 
For /c = 0, . . . , -fC let Jfc : 5fc — > M have its usual meaning; and let : rekNk — > M be the function 
defined by Ik{x) := — e|)~^''^ using the coordinates of the neck introduced in Section [3J 

This latter function is in the kernel of the linearized mean curvature operator of the catenoid with 
respect to the Euclidean background metric; it is an odd function with respect to the center of the 
catenoid and is asymptotic to ±1. Now for convenience let / := fr{(T,5) and := ij,rf{T,r{a,5)) 
denote the solution found in the previous section and define vr^ : VV* —>■ by 

^fc(e) ■■=( [ e • xZck,kIk dVol, , / e • x:tt.kJk dVol, ) . 



The notation for the cut-off functions from Section [6TL2] has been used here. A consequence of the 
following lemma is that if 7r(e) = then e = 0. The proof is a straightforward computation. 

Lemma 26. Choose e £ W and write e = Ylk=i {'^kx2xt,kJk + bkXllt,k^s{vtxZck,k)) f"^ 
Ofc, 6fc £ I^- Then 

7rfc(e) = {Cibk + C'ir'^{e^J'^ak - > C2r'^ak) 

where Ci,C[,C2 are constants independent ofr, e and 5. 

A fundamental application of the expansions of the mean curvature found in Proposition [1] 
and equation ([7|) from Section l2.ll is to derive a formula relating ■n{£{wr{(y,li)) to the geometry 
of S,.(cr, (5). It is via this formula that the location of the spheres and necks in I],.(cr, (5) and the 
background geometry of M conspire to determine when a nearby CMC surface can be found. 

Proposition 27. Let := /ij(Sr(cT, 5)) and f := fr{(T,5) for convenience. Then the mean 
curvature of T,f satisfies the formulce 



I - f ) xZcKuh dVol, = C,5ure^l^ 



(29a) 

+ 0(eii/4r2) + 0(r2e3/2) + 0{r^R{r, e, 6)e^) 
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JT,f (29b) 

+ 0{er^) + 0(rV/2) + 0{r^R{r, e, 6)) 

where Ci,C[,C2 are constants independent of r, e and 6; and S{pk) ■= 5'(V^(-(^)5'(pfc), 7(tfc)) is the 
component of the gradient of the scalar curvature of the ambient metric along the geodesic 7 at pk- 

Proof. The formula (|29bp will be derived first. Consider the subset Vk of consisting of the k*'^ 
spherical region of S := Sr('7) ^) and its adjoining transition regions. Let X = and X = Xg*j j^X 
using the normal coordinate system centered at pk used in the definition of S n V^. Furthermore, 
define the domain in M by requiring dWk = VfcUciUc2 where ci and C2 are small embedded two- 
dimensional disks with boundaries dci and dc2 contained in t = constant planes with ^r-conormal 
vectors i^i and 1^2 tangent to Vk- 

Now let T,' be either n Vfe or S n V^; let X' be a vector field supported on this surface; and 
let g' be any choice of background metric. Define I{T,' , X' , g') to be the integral in (|29bp except 
with Tif replaced by T,' and replaced by g'{X',N') where N' is the (^'-unit normal vector field 
of S' and the mean curvature and volume form calculated from g' . It is now simple to phrase the 
means by which the formula ()29bp will be found. First, I(T,f,X,g) can be expressed as 

I{^f,X,g) = I{^,X,g) + {l{^,X,g) - I{^,X,g)) ^^^^ 
+ {l{^f,X,g) - I{^,X,g)) + • 

Then one can apply the first variation formula in Euclidean space to the first term, yielding a pair 
of boundary integrals; one can apply the expansions for the mean curvature with respect to the 
perturbed background metric from Lemma [1] to the second term, yielding a curvature quantity; 
and one can treat the third and fourth terms as small errors. 

The details of the computation outlined above are as follows. For the first term, the classical 
first variation formula for a surface with boundary in Euclidean space gives 

/ (iJ[S]-f)jfcdVolo = V(-l)^f / 5(X,i>,)dLo + f / 5(X,|)dVolo') 

= j2{-iyr{Cre, + C[eT) 

s=l 

by direct computation, where is the neck scale parameter associated to the neck adjoining the 
curve dcs, while Ci and C( are constants independent of r, e and 5. For the second term, the 
expansion in equation [7] implies 

I{^,X,g)-I{Y.,X,g)= [ iV) ((iRic(y, + j^VyRic(y, F)) - f ) 



- {lRm{Ei, Y, Ej,Y) + iVyRm(Si, Y, E,,Y))B'^ 

- fRic(iV,y) - iVyRic(iV,y) 

+ ^Vj^mc(Y,Y) - iV^Rm(iV,y,7V,y))dVolo + 0(r 



(31) 
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where Y is the position vector field of Sjt, while the quantities g, N, H, B and N, H, B have their 
usual meanings. Since Vk is the normal graph of the function rG over the sphere as in Section 
[3l one can replace the integral in (j3ip with an integral over 5^, at the expense of an error of size 
0{£\ log(e)|r^). Hence by direct computation using H = ^ and Bij = rhij one finds 

/(S, X, g) - /(S, X, g) = -C2r^S{pk) + 0{e\ log{e)\r^) 

where C2 > is a constant independent of r, e and 6. An similar computation is performed in j2Uj . 

It remains to estimate the error terms appearing in (jSOp . In the third term, the fact that 
H[^f] = H[T] + C{rf) + Q{rf) must be used. Thus 

using the estimate of / from Proposition [20l In the fourth term, observe that X — X is supported 
in a collar of width 0{r£^l^) around the transition regions of S Pi V^. Hence 

|/(S,X,5)-/(S,X,<7)| <CrV/2 

using the estimate from Step 3 of Proposition [13] for the mean curvature in the transition region. 

The derivation of the formula (j29a|) is similar to what has been done above. That is, writing 
I{T,' , X' , g') as before, but with X' equal to either X := or X := x^neck kM^ finds the same 
decomposition as (150]) for I{T,f,X,g). But now, 

f {Am - l)h dVolo = -- I Ik dVolo = C4r4^' 

where 5k is the displacement parameter of the neck J\fk and C is a constant independent of r, e 
and 5. This is because /j]^^^! ^[5^] /fc dVolp = exactly (this is the first variation formula for 
the exactly minimal surface S n supp(x^gcfc A;)) function with respect to the neck 

having 5^ = 0, whereas the integral is being taken over the neck with 5^ 7^ 0. The remaining terms 
in the expansion of I{Y^f,X,g) are small error terms whose estimates are sufficiently similar to the 
analogous ones above and will not be repeated for the sake of brevity. □ 

7.3 Proof of the Main Theorem 

The formulae developed for the balancing map Br in the previous section make it possible 

to choose an exactly CMC surface from amongst the family of surfaces 5) (^^(cr, (5)) . This 

will be done as follows. First, because of Lemma [26l it is sufficient to find {cr,6) so that the right 
hand sides of equations (j29ap and ()29bp vanish for every k. One should realize that pk in these 
equations can be expressed in terms of £1,621 • • • via the formula pk '■= j{2kr + X]f=i ^i) ^^^d £k 
can be expressed in terms of ak via the formula ak ■= ^ki^k) as in Section [3l (In the case of the 
one-ended surface, let the relationship £k = Pt{0) for T = 2 + crx/r satisfied by the Delaunay end 
of T,^^{a,6) be written sk ■= Ak{ctk) for consistency.) Note that ak and are smooth functions 
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of Ek- Finding the appropriate value {(7,5) will amount to applying the implicit function theorem 
for smooth maps between finite dimensional spaces to this system of equations, and will lead to a 
unique solution {cr,5) := {a{r),6{r)) for all sufficiently small r > and e, 6 satisfying < e < 
and 5 < e^/^. 

The finite-length surface. The equations that must be solved to produce the finite-length CMC 
surface are as follows: if the various error quantities appearing in equations ([29]) . divided by r, are 
denoted Eg^k{r,e) where s = 1 refers to a neck and s = 2 refers to a sphere, then 

: (32a) 

= 6k-i + e]^^^'^Ei^K-i{r,e,6) 

as well as 

= -qiiei) + q2{e2) - C2,ir^S{pi) + E2,i{r,£,5) 

= -Qkiek) + qk+i{£k+i) - C2,kr^S{pk) + E2,k{r,e,6) (32b) 

= -qxieK) - C2,Kr^S{pK) + E2,K{r, e, 6) 

where gfc(e) := Ci^^e + C{ ^e'^/^ and pk := j{2kr + J2t=i^i) while Ci^/c,C{^ and C2^k are various 
constants independent of r, e and 6. Note that the Eg^k are smooth functions of e. Also, because 
the t —t symmetry that has been imposed since the beginning, the scalar curvature must have 
a critical point at po- 

One should now view the equations in (|32|) as a systems of equations for the e and 6 variables 
depending on the parameter r that is to be treated using the implicit function theorem. When 
r = there is an exact solution 6i = ■ ■ ■ = 5k-i = and £i = ■ ■ ■ = ek = 0- Furthermore, it is 
easy to see that the derivative matrix of the function <I>(e,(5, r) defined by the right hand sides of 
(j32ap and ()32bp in the e and 5 variables is invertible at r = with a lower bound of size 0(1) on 
its determinant (the derivative matrix is upper-triangular with non-zero constants of size 0(1) on 
the diagonal). Hence by the inverse function theorem there is a solution of (j32p for all sufficiently 
small r, and the dependence of e and 5 on r is smooth. Note that the solution for small r will 
have Ek = C){r'^ Ylk'=o ^(Pk')) and hence Cir^ < Sfc < C2r^ for numerical constants Ci and C2. 
This is because the sum r^^,^Q S{pk')) approximates a Riemann sum for the integral of S along 
7 from t = to t = 2Kr and a uniform bound on the oscillation of the scalar curvature of the 
ambient manifold has been assumed. Furthermore, it is also the case that 5k < e^^'^ for small r 
simply by examining the dependence of the E(r, e, 6) quantities on its arguments. This completes 
the construction of the finite-length CMC surface. □ 
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The one-ended surface. The equations that must be solved to produce the one-ended CMC 
surface are shghtly different. Using the same notation as above, these equations are 

: (33a) 
= 6k-i + e^^^'^Ei^K-i{r, e, 6) 

as well as 



= qo{eo) - C2,or^S{po) + E2,o{r, e, 6) 

= -qoieo) + qiiei) - C2,yS{pi) + E2,i{r,e,5) 



(33b) 



= -QK-iisk) + QKieK) - C2,Kr^S{pK) + E2,K{r, e, 5) . 
One should again view the equations in (|33p as a systems of equations for the e and 6 variables 
to be treated using the implicit function theorem, but this time depending on the parameters r and 
the point po. When r = and pQ is any point on 7, there is an exact solution 5i = ■ ■ ■ = 5k-i = 
and El = ■ ■ ■ = ek = 0- Furthermore, the derivative matrix of the function of (e, 6, r) defined by 
the right hand sides of (j33p in the (e, 6) variables is invertible at r = with a lower bound of size 
0(1) on its determinant. Hence by the inverse function theorem there is a solution of (|32p for all 
sufficiently small r, the dependence of e and on r is smooth, and the dependence of the solution 
on r is the same as before. This completes the construction of the one-ended CMC surface. □ 
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